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Abstract. We propose a simple adaptive-network model describing recent swarming 
experiments. Exploiting an analogy with human decision making, we capture the dynamics of 
the model by a low-dimensional system of equations permitting analytical investigation. We 
find that the model reproduces several characteristic features of swarms, including spontaneous 
symmetry breaking, noise- and density-driven order-disorder transitions that can be of first or 
second order, and intermittency. Reproducing these experimental observations using a non- 
spatial model suggests that spatial geometry may have a lesser impact on collective motion 
than previously thought. 
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1. Introduction 

"More is different," a central observation in complex systems research is in few areas as 
evident as in collective decision making. Recent studies on groups of self-propelled agents, 
such as bird flocks, fish schools, insect swarms, or herds of quadrupeds show that these can 
often make better choices in groups than individually [1-3]. Effective collective action, which 
has given rise to the popular notion of swarm intelligence, appears to follow from universal 
organizing principles [4]. The ease and elegance with which, e.g., a school of fish finds its 
way and avoids predators, can be even more appreciated if one compares it for instance to 
political decision making in humans. 

It is interesting to note that, although regarded as collective decision-making processes, 
swarming behaviour and collective motion are modelled differently from decision making in 
human populations. Theoretical studies of collective motion have mostly focused on agent- 
based simulations of self-propelled particles [5-10] and fluid-like models that treat them as 
continuous media [11, 12]. By contrast, studies of decision making and opinion formation in 
social systems typically represent the system as a network, emphasizing the discrete nature 
of interactions [13-15]. One of the reasons for this difference in modelling approaches is 
that for swarm systems, spatial embedding is assumed to be of central importance, whereas 
social interactions are felt to be less constrained by physical space. Furthermore, in collective 
motion, it cannot be neglected that an agent's decision to move in a certain direction 
determines the agents with whom it will be interacting next. However, a similar feedback of 
individual decisions on future interaction partners was also studied in recent works on opinion 
formation [16-23]. The resulting models incorporate both an opinion formation process on 
the network and a dynamic update of the network topology and thus fall into the class of 
adaptive networks [24, 25]. 

In the present article, we propose a non-spatial adaptive-network (AN) model of 
swarming behaviour. Reproducing characteristic observations for swarm systems, we find 
that spatial geometry might play a less central role than has been assumed. Our approach 
highlights the analogy between swarming and social consensus, thus building a bridge 
between two areas of research that have been so far considered separately. 

2. Adaptive-network model for a swarming experiment 

We focus on the swarming experiments by Buhl et al. [26]. In their set-up, groups of 5 to 
120 locusts were placed in a ring-shaped arena and left to march freely for 8 hours while a 
digital camera captured their positions and orientations. At low insect number, no ordered 
collective motion arises; the system displays no clear clockwise or counter-clockwise flow of 
locusts around the arena. At intermediate insect numbers, locusts start aligning, generating 
long periods of collective rotational motion during which most agents are marching in the 
same direction. These periods of coherent motion are interrupted by rapid spontaneous 
changes in their collective heading direction. Finally, at high insect number this spontaneous 
direction switching is no longer observed and agents rapidly adopt a common and persistent 
marching direction, either clockwise or counter-clockwise. Buhl et al. reproduced these 
experimental results qualitatively in simulations using a one-dimensional agent-based model 
of self-propelled particles, also investigating the effect of inherent noise in a subsequent study 
[26, 27]. 

Here, we model this experiment using a different approach, trying to address the 
mechanism leading to the observed collective dynamics with the help of a simple low- 
dimensional description that lends itself to analytical treatment. To that end, we consider 
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the system of interacting agents as a complex network. Each node represents an insect, and 
nodes are linked if the corresponding locusts are mutually aware of each other through any 
interaction mechanism. As in [26], we distinguish only two directions of motion: every node 
can be in an R or L state, representing an agent that marches clockwise (a right-goer) or 
counter-clockwise (a left-goer), respectively. We refer below to pairs of nodes (agents) in the 
same state as equal-goers and to those in different states as opposite-goers. 

The proposed model only takes into account the agents' headings and contact network 
while neglecting all other information, including insect positions. The evolution of the 
network is modelled by a set of stochastic processes. As locusts advance in the experimental 
system, non-interacting opposite-goers eventually meet and start sensing each other. We 
model this by randomly introducing R-L links at a rate of a per node. In addition, interacting 
opposite-goers will eventually lose contact, which is modelled by the random deletion of 
R-L links at a rate of d a per link. Likewise, equal-goers can start or stop interacting as 
they approach or separate from each other due to marching speed differences or lateral 
displacements. This is represented by also introducing for equal-goers the attachment rate per 
node a e and deletion rate per link d e . Using these conventions, all rates are defined as intensive 
quantities. The state dynamics of each node is given by a stochastic process that depends on its 
topological neighbours. We assume that each node switches direction with probability W2 for 
every R-L link it has to an opposite-goer. To account for non-linear three-agent interactions, 
we introduce an additional probability W3 of the central node switching direction for every 
L-R-L and R-L-R chain. Finally, noise is represented by a constant probability q of an agent 
spontaneously switching direction. 



3. Analytical description 

In order to study the collective dynamics of the AN model, we define convenient observables, 
the so-called moments, given by the densities of various subgraphs in the network [28, 29]. 
Each subgraph can be classified by its order, i.e. the number of links it contains. Zeroth 
order moments are given by the right- and left-goer densities ([R] and [L], respectively). First 
order moments are the per-capita densities of R-R, R-L, and L-L links ([RR], [RL], and 
[LL]). Second order moments correspond to the densities of A-B-C triplets [ASC], with 
A, B, C £ {R, L}. The moment dynamics is captured by balance equations containing the 
variables of interest together with densities of larger subgraphs. 
The zeroth-moment equations are 

±[R] =q([L]-[R}) + w 3 ( [RLR] - [LRL]) , (1) 

and the symmetric expression for [L], obtained by interchanging R and L. The first-moment 
equations are 

^[RR] = q{[LR] - 2[RR]) + w 2 {[LR] + 2[RLR] - [RRL]) 

+ w 3 (2[RLR] + 3[ R Lft] - [ R R L L ]) + a e [R} 2 - d e [RR], (2) 

and the symmetric expression for [LL\. Here we use [^B^] to denote the density of 
third-order motifs with a central node in state B £ {L, R} linked to three nodes in states 
A, C, D £ {L, R}. Finally, rather than writing an equation for [LR], we note that the total 
first-moment link dynamics depends only on the link creation and deletion processes through 

— ([LR] + [RR] + [LL]) = a [L][R] - d [LR] + 

a e ([R] 2 + [L] 2 ) - d e ([LL] + [RR]) . (3) 
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Figure 1. Top: Bifurcation diagram of the density of right-goers [R] vs. link creation rate 
a . Solutions of the ODE system of Eqs. (l)-(3) (solid line) yield a supercritical pitchfork 
bifurcation in excellent agreement with results from numerical network simulations (circles). 
Bottom: Phase diagram showing the bifurcation point as a function of the link creation 
rates a and a e . In the bistable region (grey), the pitchfork bifurcation becomes subcritical. 
Parameters: N = 10 4 nodes, d = 0.25, d e = 0.1, 102 = ^3 = 0.2, q = 0.1, and (top 
only) a e = 0.2. 



The ODE system (l)-(3) can now be closed using a pair-approximation [28-30], where 
triplets and quadruplets are given by 

[RLR] = J-^-, [i^H2.Mra, (4) 



rRrRl .,2 [ LR f rflpil ,„2 [LR] 2 [RR] 



and symmetric expressions. The factor k = ((k 2 ) — (k)) / (k) 2 relates the second and first 
moments of the degree distribution. Because our network dynamics will yield an unknown, 
randomly evolving topology, we use a random graph approximation, setting k = 1 as in 
[18,21,28-30]. 

For a e = d e = 0, the stationary solutions of this ODE system decouple, with Eqs. 
(1) and (3) solved independently. We obtain analytically a mixed-phase solution branch 
([R] = [L] = 1/2) that becomes unstable in a supercritical pitchfork bifurcation at 
a* = 2d y / 2q/ (kw^), giving rise to the two collective-motion solution branches [R]± = 
(1/2) ± y/1 — Sqd^/ (nw^a 2 ) /2. For a e ^ and d e ^ 0, the stationary solutions can be 
computed numerically by solving the corresponding system of algebraic equations. Here, we 
also find a supercritical pitchfork bifurcation for small a e , as shown on figure 1 (top). But for 
higher values of a e , the transition occurs through a subcritical pitchfork bifurcation (figure 2, 
right-column insets). This yields a bistable phase where ordered and disordered states coexist, 
highlighted in figure 1 (bottom). 

We note that a supercritical pitchfork bifurcation, given by [R]± = (1/2) ± 
i/l — Aq/wz/2, could be already observed in a simpler approximation, in which the system is 




Figure 2. Degree distributions of four different stationary solutions obtained through adaptive 
network simulations. The top (bottom) row shows cases in the disordered (ordered) phase, 
with insets displaying their location in the bifurcation diagram. The left (right) column shows 
cases with a supercritical (subcritical) pitchfork bifurcation. Different curves display the 
connectivity of left-goers (L), of right-goers (R), of left-goers only to left-goers (LL) or only 
to right-goers (LR), etc. Parameters: Same as in figure 1, except when noted on plots. 



closed at the zeroth order. However, the pair approximation is more accurate when compared 
with stochastic simulations of the network. Moreover, it allows for the new class of subcritical 
solutions, which we discuss below. 

4. Results 

In figure 1 (top), we show that the ODE system solutions are in excellent agreement with 
individual -based stochastic simulations of the network dynamics. We verified that the small 
remaining discrepancy is due to the k = 1 assumption, and not other factors such as finite-size 
effects. 

We now compare the results of our AN model to the locust experiments in [26]. The 
bifurcation diagram in figure 1 uses the encounter rate between opposite-goers a as control 
parameter, which is proportional to the experimental agent density in the current framework. 
We can thus compare our results directly to the behaviour observed experimentally at different 
agent numbers. For a low association rate a , the network exhibits no order or symmetry 
breaking. This corresponds to the disordered mix of left- and right-goers showing no 
collective motion obtained in experiments at low insect number. 

For high a a , the system must be in one of the two solution branches, with a majority of 
nodes in either state (R or L). This is the ordered collective marching state found at high insect 
number [26]. A similar pitchfork bifurcation is also observed when using the noise intensity, 
q, instead of a as control parameter, as done in most previous numerical works [5, 6, 9]. 

Let us emphasize that the presence of this transition in the AN model implies that very 
few elements of the agent dynamics are required to obtain such swarming behaviour. In 
particular, we did not choose any specific interaction rule but only required that it drives 
agents to head in the same direction. By contrast, we find that three-body interaction processes 
are required to break the symmetry and obtain swarming solutions. Furthermore, a subcritical 




Figure 3. Cumulative distribution of residence-times in a majority L or majority R state 
for the adaptive network dynamics at various values of a a . As a is increased, the system 
becomes more ordered, switching direction less often, and the distribution becomes broader. 
For a > 0.6, a preferred residence time appears at large r ~ 10' 5 due to finite-size effects. 
Insets show the density of right-goers [R] vs. time for a a = 0.1 (bottom-left) and a = 1.0 
(top-centre). Parameters: Same as in figure 1 (top) but for 100 nodes. 



bifurcation, giving rise to hysteresis or sudden polarisation, is only possible in the AN model if 
a e and d e are non-zero. This qualitative result could shed light on the current controversy over 
the order of the swarming transition [5, 6, 31-33]. Indeed, a first-order transition (stemming 
from the subcritical pitchfork bifurcation) is only possible here if groups of equal-goers can 
associate or dissociate while heading in the same direction. It would be very interesting to 
explore if a similar effect is present in agent-based simulations and experiments. 

Figure 2 shows the degree distributions obtained in the ordered and disordered phases 
for supercritical and subcritical cases. Both display similar connectivities. In the disordered 
states (top), most agents have very few links to equal-goers and no links to opposite-goers. 
This is also observed in agent-based simulations and experiments [9, 26], where the disordered 
regime develops no large clusters and, therefore, small connectivity. 

In the ordered (right-going) state (bottom), R-R links are strongly favoured. This 
corresponds to the formation of large right-going groups in the agent-based dynamics. The 
number of L-R links also increases, which corresponds to encounters between a few left-goers 
and these large right-going clusters. The typical number of all other links decreases. 

In previous experiments and simulations it was observed that there is an intermittent 
regime where the swarm is polarised but can switch the marching direction spontaneously 
[26, 34]. This behaviour is also captured by the AN model. The intermittent regime occurs at 
intermediate densities, i.e. at a a values close to the bifurcation, in the ordered phase, where a 
low nucleation barrier allows for stochastic switching between the two branches. 

Figure 3 shows the cumulative distributions of residence times (lasting r or longer) 
in which the network resides in a majority R state before switching to a majority L state 
or vice-versa. In the disordered phase, this distribution decays exponentially as expected, 
since it results from memory-less stochastic fluctuations about the stationary state. As a 
approaches its critical value, the distribution develops a long tail approximating a power law 
with exponent —1/2, providing evidence for a switching process with memory. We note 
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that a power-law distribution of switching times with the same exponent was observed in 
previous spatial models [9, 34]. The experimental locust dynamics also display intermittency 
at intermediate densities, but the available time series are not long enough to characterize its 
statistics [26]. 

For a > 0.8 finite-size effects produce a preferred residence time at large r ~ 10 3 
that grows with the system size, appearing as a plateau in the cumulative distributions. This 
corresponds to the typical escape time from highly polarised states where the system gets 
trapped when system-wide connectivity is reached within the finite network. 

5. Conclusions 

In summary, we have proposed an adaptive-network model of a swarm experiment that 
captures much of its characteristic collective behaviour and likens it to an opinion-formation 
process. In particular, our model displays a transition from a disordered to an ordered phase 
with increasing insect density. Furthermore, an intermittent regime is observed close to the 
transition point, where a fat-tailed distribution of residence-times emerges. 

We emphasize that in the proposed model, these characteristics of swarming systems 
are recovered without an explicit spatial representation of the system. This suggests that 
the spatial context of swarming may not be of central importance for many phenomena. By 
contrast, three-body processes and an increased probability of interaction between agents with 
intersecting trajectories are found to be essential in our model. 

Further work is certainly necessary to test whether the experimentally observed transition 
is caused by the same mechanism that is at work in our non-spatial model. In order to improve 
our understanding of the role of spatial dynamics in swarms, detailed network- and agent- 
based simulations should therefore be compared. We postulate that the same dynamics will 
be followed by point-like agents, which are well described by the nodes in our model, but not 
by spatially extended agents, where jamming must play an important role. 

In the present work, we used a modelling approach originally proposed for social 
networks. We believe that the analogies we have drawn between swarming phenomena and 
opinion formation processes could be fruitfully exploited in further studies, thus building a 
bridge between the two fields. 
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